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Abstract: We exposit the eigenvalue distribution of the lattice Dirac operator in Quan- 
tum Chromodynamics with two colors (i.e. two-color QCD). We explicitly calculate all the 
eigenvalues in the presence of finite quark chemical potential \i for a given gauge configu- 
ration on the finite-volume lattice. First, we elaborate the Banks-Casher relations in the 
complex plane extended for the diquark condensate as well as the chiral condensate to 
relate the eigenvalue spectral density to the physical observable. Next, we evaluate the 
condensates and clarify the characteristic spectral change corresponding to the phase tran- 
sition. Assuming the strong coupling limit, we exhibit the numerical results for a random 
gauge configuration in two-color QCD implemented by the staggered fermion formalism 
and confirm that our results agree well with the known estimate quantitatively. We then 
exploit our method in the case of the Wilson fermion formalism with two flavors. Also we 
elucidate the possibility of the Aoki (parity-flavor broken) phase and conclude from the 
point of view of the spectral density that the artificial pion condensation is not induced by 
the density effect in strong-coupling two-color QCD. 
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1. Introduction 



Quantum Chromodynamics with two colors (two-color QCD) instead of three is a sophis- 
ticated practice ground for theorists to extract worthwhile information out of dense quark 
matter. We immediately hit on several reasons why we can believe so: First of all, nu- 
merous works on dense two-color QCD have almost established a firm understanding on 
the ground state of two-color QCD by the analytical approach as well as the Monte-Carlo 
simulation § |, § | | g @, §, | |HJ 0, 0, |1, 0, 0, M, 0, B © 11,11,11, 
25, 0, |28], HH, |3^, ^1J| . Second, the notorious sign problem of the Dirac determinant 



at /x ^ (where is the quark chemical potential) is not so harmful as genuine QCD, 
which makes it viable to perform the Monte-Carlo integration [|], [ll], 0, |3^| . Third, dense 
two-color matter realizes a bosonic baryon system leading to the Bose-Einstein condensa- 
tion of the color-singlet diquark ||, fjj. This two-color superfluid phase is reminiscent of 
the three-color superconducting phase p3[ ], for they both break the Ub(1) symmetry. Fi- 
nally, enlarged flavor symmetry earned by the pseudo-real nature of the SU(2) group that 
is called Pauli-Giirsey symmetry p5| , 34] constrains two-color QCD at m q = fj, = 0. The 
interplay between the chiral and diquark sectors simplifies owing to the symmetry, which 



enables us to construct an effective model for two-color QCD with less ambiguity 21, 25]. 

This paper aims to illustrate the spectral behavior in a rather brute-forth manner. 
We usually define the order parameter and concern its expectation value to examine the 
phase structure with varying the external parameters such as the temperature T, the quark 
chemical potential fi, the quark mass m q , and so on. We shall explore our another trail 
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here leading to the phase distinction. In this work we will carefully look into the eigenvalue 
distribution of the Dirac operator and characterize the state of matter by the distribution 
pattern. Actually, we can find the scatter plot for the Dirac eigenvalues in two-color QCD 



in Refs. 0, 11] and we will do that in a more systematic way. It is long known that the 



eigenvalue spectrum is informative in the vacuum 36] and the random matrix theory is 
capable of determining the low-lying spectrum, which has recently been extended to the 
finite density study 31, 37|. Interestingly, the comparison to the random matrix model 
exhibits good agreement also in the case of the overlap fermion at fi ^ pq] . 

It is not only the low-lying spectrum but also the whole spectral density that we will 
deal with in the present paper. The Monte-Carlo simulation generates a set of gauge 
field configurations each of which has a substantial weight on the partition function. One 
configuration corresponds to one value for a certain operator (the order parameter for 
example of our interest), and the more configurations we accumulate, the more accurately 
we can improve the expectation value of the order parameter. Here, we would remind that 
the well-known Banks-Casher relation |39| yields the chiral condensate given in terms of 
the eigenvalue spectral density at the origin (i.e. m q — ► 0). It follows in turn that the 
order parameter makes use of only tiny amount of the entire information available from 
the spectrum. In this work, hence, we will unveil detailed information in a special case of 
dense and cold (T = 0) quark matter with two colors. 

One might come across a question then; what is the benefit from the whole shape of 
the eigenvalue distribution? To answer this, we should be aware that the Dirac eigenvalues 
originally lie on the imaginary axis except for the displacement in the real direction by 
m q but they scatter over the complex plane because of non-zero chemical potential fi ^ 
or Wilson coefficient r ^ 0. This feature has, more or less, something to do with the 
sign problem meaning that the Dirac determinant could take a negative value. For [i ^ 
the Dirac operator mixes the Hermitean and anti-Hermitean operators up resulting in a 
complex eigenvalue. The situation at finite \x looks similar to that in the presence of the 
Wilson term in view of the eigenvalues particularly in the two-color case |3^] . [We implicitly 
assume only the two-color case below.] The sign problem may arise actually when the 
eigenvalue distribution protrudes from the positive quadrant into the negative quadrant. 
This observation implies that a large m q (center location of the eigenvalue distribution) as 
compared to the chemical potential or the Wilson term (distribution width) would put the 
sign problem away. In the physics language, the vacuum stays empty as long as m q > fx, 
and so there is no \x dependence then, which brings about no sign problem naturally. The 
onset of the density effect is manifestly visible from the whole eigenvalue distribution. 
Besides, since m q ^ shifts the distribution, it is transparent to take account of the mass 
effect provided that the eigenvalue distribution is given. These motivate us to turn to the 
entire Dirac spectrum. In the future, hopefully, we believe that the eigenvalue distribution 
should shed light upon the sign problem at a deeper level. In fact, as we will recognize 
later, a large value of fi induces a peculiar structure in the eigenvalue pattern. 

The above mentioned may well be somewhat abstract. Let us then make the issue to 
be discussed more specific. What puzzles us is that there seems to be no clear distinction 
between the onset criterion for the superfluid phase and the Aoki phase if considered based 
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on the eigenvalue distribution alone. They can possibly coexist but it would be a weird 
situation because the superfluid phase is a physical ground state but the Aoki phase is a 



lattice artifact inherent in the Wilson fermion formalism [40, 41]. The final part of this 
paper will be devoted to resolving this matter. There, we will find that the onset criterion 
is certainly degenerate when [i = 0, but the Aoki phase is taken over by the superfluid 
state in the proper limit of fj, — > 0. In short, we conclude that the Aoki phase never 
emerges by the density effect in strong-coupling two-color QCD. This statement does not 
conflict the preceding strong coupling analysis |4(], ^l| because the Aoki phase solution at 
strong coupling is a saddle-point and infinite N c is required for stability, though this fact 
is sometimes overlooked. 

2. Two-Color QCD at Strong Coupling 

In the limit of the strong coupling the gauge action does not enter the dynamics and the 
partition function is simply given by the fermionic part; 

(detV) Nt } = fY[dU„(n)(detV) N! . (2.1) 

Here T> is the Dirac operator. Although the strong coupling limit is a drastic approxima- 
tion which neglects the gauge dynamics completely, it is amazing that only the Dirac de- 
terminant with random gluon fields can grasp rich contents of quark matter not only in the 
two-color case @, §, |§ but also in the general case g g(3[ gl|, g|, @ @> @ © 0, EH EH" 



In a box with volume V = L 4 , the operator P is a (4^^) x (4N C V) matrix. We 
denote the eigenvalue of V by Aj, that is, 

Vvi = X iVi , (2.2) 

where i runs from 1 to AN C V . Then the Dirac determinant is given by the product of all 
the eigenvalues. It is easy to prove that detP in the SU(2) gauge theory takes a real value 
even at finite density where V loses the 75-Hermiticity, i.e. 7515(^)75 = T>>(—[i) ^ V^(n). 
The standard argument immediately follows; 

det £>(//) = det[(Ca 2 j 5 y 1 V(fi)(Ca 2 ^)] = detPV) = [det £>(//)]* . (2.3) 

Here, to derive the above, the necessary relations we use are 757^75 = —7^, Cj^C^ 1 = 
— 7J, and = U* where the last relation corresponds to the pseudo-real nature of the 

SU(2) group. 

From this argument we see that detD(/u) is real but not necessarily positive. The 
simulation thus entails an even number of iVf so that (detP)^ is positive definite. This 



is the main reason why the exotic phase structure proposed in Refs. [12, 25] in two-color 
QCD with quark and isospin chemical potentials has been far from confirmed. The two- 
color determinant, however, buries a nice property of respective eigenvalues under the 
product. We can prove that, if Aj = m q + iA^ is an eigenvalue of the Dirac determinant 
in two-color QCD, there appear m q — iA^, m q + and m q — simultaneously in the 
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eigenvalue spectrum [11, 25, |32]. The proof may break down when iA- is a real number; the 



eigenvectors for m q + iA^ and m q — \\'* could not be independent. According to Ref. [11] 
the staggered fermion is safe from such a possibility but the Wilson fermion has only a 
pair of m q + iA^ and m q — iA^ instead of a complex quartet in that case of real iA ■ . We will 
explicitly verify that this is the case. Then the single-flavor Wilson fermion suffers the sign 
problem once either of real m q + iA^ and m q — iA^ is negative. 

Here we shall briefly summarize the known facts in two-color QCD at strong coupling. 
Let us begin with the chiral limit. It has been discovered first in Ref. Q that the chiral 
condensate is zero, while the diquark condensate has a finite expectation value, in the 
limit of jj, — > with m q = taken first. In the presence of m q ^ the system is kept 
intact as long as fi is sufficiently small and, in turn, the chiral condensate becomes non- 
zero but the diquark condensate vanishes. As soon as \x exceeds the mass of the lightest 
excitation (usually bosonic baryon), the density effect is activated leading to decreasing 
chiral condensate and increasing diquark condensate as fj, goes larger. We remark that 
this behavior of dense two-color QCD has been settled in the staggered fermion but the 
relation between the chiral and diquark condensates is not quite convincing yet in the 
Wilson fermion because the Wilson term breaks chiral symmetry explicitly. 

3. Banks-Casher Relations 

Here we will make a quick view over the link between the eigenvalue spectral density and 
the chiral, diquark, and parity-flavor breaking condensates for later usage. In this section 
the argument holds regardless of strong coupling or not. 

3.1 Chiral Condensate 

It is widely known that the chiral condensate has a close connection to the Dirac eigenvalue 



distribution via the Banks-Casher relation [39|. To advance our discussions in a self- 



contained manner we shall take a brief look at the derivation of the Banks-Casher relation. 
In the explicit presence of the source for the chiral condensate (i.e. mass term), the Dirac 
operator could be decomposed into the form of T>[m] = m q \ + T>[0] whose eigenvalue is 
denoted as A» = m q + iA^ as we did in the previous section. The chiral condensate per 
flavor is given by the derivative of Z with respect to m q , which leads us to 



(3.1) 



where p x (X) is the eigenvalue spectral density which is to be expressed in the complex 
plane as 

"*M = ^Evb;> <M) 

i 

which is, strictly speaking, the resolvent [^] rather than the spectral density. To keep the 
analogy to the conventional Banks-Casher relation, however, we shall refer to the above 
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as the spectral density. The integration contour should go around all of the poles at Aj to 
pick all the eigenvalues up. In our notation (•••){/ means the ensemble average over gauge 
configurations and ((•••)) represents the average including the Dirac determinant. 

Here we consider the contour which is an infinitely large circle in the complex plane 
surrounding all the poles. Then the contour integral must amount to zero because p x (X)/\ 
goes to zero faster than |A| _1 . That means that we can evaluate the above integral by the 
negative residue of the pole at A = 0. After all, we have 

<^> = -N f 7T«p x (0))) . (3.3) 

For consistency check let us consider a bit more about this formula. Usually the Banks- 
Casher relation is given in the limit of m q — > 0. In this limit, using the notation \ = m g +iA^ 
where A^ is a real number in the continuum theory, we can rewrite Eq. (|3.2|) into a form of 



p x (0) = lim \ Y, — = \? E b ^ > ( 3 ' 4 ) 
m.—»o ttV m + iA- V 



which is more familiar in literatures. We note that Eq. ( |3.2| ) is in fact an analytic continued 
form of the expression Q3.4| ) with the delta function, and it is equivalent to the definition of 
the resolvent used in the context of the random matrix theory |2J| . This complex extension 
is necessary for our purpose since the Dirac operator loses Hermiticity at finite density or 
in the Wilson fermion formalism. One might have noticed that the Banks-Casher relation 
in the complex plane is a trivial relation; it is obvious from Eq. (|3.2| ) that — irp x (0) returns 
to Yli \ rl immediately. 

3.2 Diquark Condensate 

We can develop the same argument for the diquark condensate as well as the chiral con- 
densate. We shall limit our discussions to the case with degenerate two-flavor (u and d) 
quarks, and then we do not need to introduce the Nambu-Gor'kov basis. In the presence 
of the same quark chemical potential [i for u and d quarks, we can write the Lagrangian 
density down as |23| 

C = tjj u V(fi)ifj u + tp d T>(fj,)tp d - J4> u (C-f 5 )<j 2 ipd + ji>d( c 75)v2ipu , (3.5) 

where J and J are the source for the diquark and anti-diquark which are anti-symmetric 
in spin, color, and flavor. By means of a variable change by 

<Pd = iPd Ca 2 > (Pd = Co 2 $ r d , (3.6) 
it is possible to compactify the above into 

ft' ."(-:>)(:;)• 

The integration over the quark fields is then straightforward and the resultant partition 
function is given as the determinant as follows; 

Z(J) = (det faf> V ~ J \ \ = (det^V) + |J| 2 ] ) u , (3.8) 
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where we have used -f^D{— y)^ = V'(fi). We note that V{y)V\y) is always Hermitean 
though V{n) may not be so. We can then prove that the eigenvalue of U(fJ,)T>> is non- 
negative real, which we denote by £f with choosing £j > 0. The diquark condensate thus 
reads 

(MC l5 )a 2 ^) = ^jZ(J)\ j=o = ^((E^TfTp)) = *«M0)» , (3-9) 
where we have defined the diquark spectral density, 

/^(0 = ^E^-^). (3-10) 



I 



in a familiar form. It should be mentioned that we do not have to perform the analytic 
continuation this time because £j sits on the real axis. 

3.3 Parity-Flavor Breaking Condensate 

In the same way we can discuss the parity-flavor breaking condensate whose non-zero 
expectation value characterizes the Aoki phase in the Wilson fermion formalism. For two- 
flavor quarks the source term for the condensate {ip^Tstp) enters the Lagrangian as 

C = $ u V(iJ,)if) u + 7/} d V(fx)tp d + H^yry^u - ^chb^d) , (3-11) 

from which the partition function reads 

Z W = (de, + - p(Mh ° ( det [^(-,) + H>] ) v . (3,2) 

It is interesting to note that Eq. (3.12) above is reduced to Eq. (|3.8|) when fi = 0. As 



a result the parity-flavor breaking condensate seems to be degenerate with the diquark 
condensate in the absence of chemical potential. Once the finite density is switched on, 
D([i)D*(—fj,) is no longer a Hermitean operator, and its eigenvalue distribution spreads 
over the complex plane. Thus, if we define the spectral density by 

Ph(v) = -\ E (3-!3) 

with the complex eigenvalue rjf of the operator Z?(//)P' (— //) with choosing Re(r/j) > 0, we 
can write the parity-flavor breaking condensate as 

<^i75^ - fes^d) = -m((p H (iH) - p H (-W))) = 27rIm«p H (i^)» • (3.14) 
4. Eigenvalue Distribution for a Random Configuration 

In this work we will take only one random configuration as a representative instead of 
calculating the ensemble average over many random configurations. Actually the eigenvalue 
distribution for one typical gauge configuration turns out to be quite informative in our case. 
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This simplification is legitimate because each random configuration equally contributes to 
a physical quantity in the strong coupling limit. So to speak, the strong coupling theory 
is democratic and any configuration is eligible for a representative. If we are interested in 
the weak coupling regime, we would have to take an appropriate ensemble average. 

We will first proceed to the calculation in the staggered fermion formalism and make 
sure that our results agree well with known results in the mean-field approximation at 
strong coupling. After that we will adopt the Wilson fermion formalism and look further 
into the possibility of the Aoki phase. 



4.1 Staggered Fermion 

The Dirac operator at finite density in the staggered fermion formalism is 
V s (n) = m q 5 m)n + - ^2 Vi( m ) u i( m ) $m+i,n ~ U}(n) 5 m , n +i 



+ 



+ m 



(m) [e" U 4 (m) i 



m+4,n 



-m rrt 



(4.1) 



where rj^n) = (— l)"i+ n 2H hn M _i anc j criem i ca i potential is introduced as formulated 
in Ref. M- 
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Figure 1: Eigenvalue distribution for a random gauge configuration in the staggered fermion 
formalism at m q = on the 6 lattice. The distribution at m q ^ is given by a shift along the real 
axis by m q . 



The zero-density Dirac operator in the staggered fermion formalism is anti-Hermitean 
except for the mass term, so that all the eigenvalues reside on the vertical line whose 
real part is m q (see the lower- left scatter plot in Fig. [l]). The chemical potential breaks 
Hermiticity and the eigenvalue distribution has a width along the real axis as [i goes larger 
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as shown in the scatter plots for = 0.3, 0.6, and 0.9 in Fig. |T[ These figures are reminiscent 
of the scatter plots in Refs. 0|. To draw Fig. |l| we have generated a random gauge 
configuration on the lattice with a volume of V = 6 4 . Because the staggered fermion has 
only color indices, the total number of dots in Fig. [I] is 2 x 6 4 = 2592 for each plot. We 
have made use of LAPACK to compute 2592 eigenvalues numerically. 

The broadened width in the real direction has a definite physical meaning. In the 
case of m q ^ the distribution has to be shifted by m q and then the entire eigenvalue 
distribution can be placed in the positive quadrant as long as \x is small as compared 
to m q . It is hence a natural anticipation that the superfluidity has an onset when the 
eigenvalue distribution becomes as wide as it reaches the origin. This is actually the case. 




Figure 2: Real part of the spectral density, Re(/o x (A)), in the complex plane for various values of 
the chemical potential. 

From the obtained eigenvalues we can explicitly calculate the spectral density (^) to 
evaluate the chiral condensate through the Banks-Casher relation in Eq. ( |3.3j ). Because of 
the quartet pattern of the eigenvalue distribution the imaginary part of p x W 1S vanishing 
on the real axis. The chiral condensate inferred from Eq. ([^) is thus insensitive to the 
imaginary part but determined solely by the real part of the spectral density taking a real 
value. We show the real part of the spectral density (^T^) in Fig. || for various \i. It is 
remarkable that the spectral density for a random configuration looks such smooth even 
without taking an ensemble average. 

As we have mentioned, the eigenvalues and thus the spectral density with a finite m q 
can be deduced simply by a shift along the real axis by m q . Therefore, /5 X (0) appearing in 
Eq. (|3.3j ) can be read from Fig. [2] by the value at (ReA, ImA) = (— m q , 0). 

When /i = fla sharp perpendicular wall stands at Re(A) = which is responsible for 
a non- vanishing chiral condensate in the limit of m q — > while keeping [x = 0. The wall is 
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smoothened by the effect of p ^ and it is no longer vertically upright at finite density, 
which leads to an interesting observation. In fact, it is not hard to conceive from Fig. || 
that the chiral condensate becomes zero in the chiral limit at infinitesimal but nonzero p. 
This is absolutely consistent with Ref. [g]. 

We shall next evaluate the diquark condensate using the Banks-Casher relation (|3.9|). 
We will start with the chiral limit (m q = 0) and then go into the finite mass case that 
we choose m q = 0.2 here in this work. For convenience we define the integrated diquark 
spectral number, 



MO 



'd£'MO 



(4.2) 



whose slope at £ = gives the spectral density Pd(£, = 0) which is proportional to the di- 
quark condensate. Although the staggered fermion Lagrangian does not involve the Dirac 
spinor, it is not difficult to make use of the Nambu-Gor'kov representation to express the 
diquark condensate by the diquark spectral density. Since the derivation is only straight- 
forward, we will not reiterate it but skip detailed arithmetics. To summarize the resultant 
relations, we can prove that 



1 



7T 



a 



(xx) = ttP(O) 



A 



1 



7T 



ix^2X) = tPd(O) 



(4.3) 



where the extra 1/2 factor in the diquark relation comes from the square-root prescription 
necessary to cancel the doubled Nambu-Gor'kov basis. In the above we have chosen the 



same normalization as Ref. [201. 
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Figure 3: Left) Histogram of no(£) whose slope gives the spectral density pd(Q- Right) Diquark 
condensate as a function of p at m q = 0. 

It is intriguing to evaluate «£)(£) by the explicit numerical calculation for the eigen- 
values in Fig. ^ from which we can get pd(£,)- Figure || shows our results in the chiral 
limit. In this case only the diquark condensate is a non- vanishing quantity Q. We plot 
the diquark condensate in the right of Fig. || without indicating the error bar. We did 
so because, though the fitting error is small, the systematic error is large. If we change 
the working procedure to measure the slope from the histogram in the left of Fig. ||, the 
resultant diquark condensate would change too. For clarity of our numerical procedure we 
explain how we compute the slope of at the origin. We assume a functional form 
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n D(0 = a£ + 6£ 2 within the range £ G [0, 0.1] and fix a and b to fit the data. Then, a gives 
the slope at the origin. If a turns negative, that means no spectral density at the origin, 
and so the diquark condensate should be zero. In this way we draw the right of Fig. || 
which shows outstanding agreement with the upper-left of Fig. 1 in Ref. j2(| . 

The m q dependence in X>(/^)X>> (//) is not such trivial as in the case of T>{fj). Roughly 
speaking, a finite m q shifts the eigenvalue in the positive real direction so that the eigenvalue 
distribution is blocked in the vicinity of the origin as long as /x is small. For /j, above a 
certain threshold value the diquark spectral density becomes finite at £ = 0, and the diquark 
condensation is activated. We can repeat the calculation in the massive case as well. As 
we mentioned our choice is m q = 0.2, and we read the chiral and diquark condensate from 
the chiral and diquark spectral density. Our final results are presented below in Fig. |j. We 
note that the onset for the chiral condensate decrease is determined by the front edge of 
the sidling wall which corresponds to the edge of the Dirac eigenvalue distribution, which 
in turn corresponds to the diquark onset. 









(1=0.7 / j-r 




j/ 11=0.95 rr" J ^ 




j-^ n=0.3 



0.05 0.1 

Eigenvalue ^ 




0.2 0.4 0.6 
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Figure 4: Left) Histogram of nr>(£) in the case of m q = 0.2. Right) Chiral and diquark condensates 
as a function of \i at m q = 0.2. 

It is impressing that the results in the right of Fig. || is consistent qualitatively with the 
mean-field analysis in the strong coupling limit given in the upper-left of Fig. 1 in Ref. p(| . 
though the direct comparison is not possible for different mass choice. 



4.2 Wilson Fermion 

We shall consider the Wilson fermion henceforth. The Dirac operator is defined as 
V w (n) = <5 m>n - «^[(r - 7,) Ui(m) 5 m+i<n + (r + 7i) u}(n)6 m>n ^ 

i 

— K (r — 74) e' 1 £7 4 (m) 6 m+ijn + (r + 74) e _/ " u\(n) 5 m)n+i 



(4.4) 



where k is the hopping parameter and we choose r = 1 throughout this work. In this 
case we adopt V = 4 4 and then there are (2 colors) x (4 spinors)x4 4 = 2048 eigenvalues. 
Of course, we could treat V = 6 4 without difficulty, but there are then many eigenvalues 
(almost five times more than the V = 4 4 case) and plotting looks too dense. Our small 
lattice volume is limited not for technical reason but for presentation convenience. 
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It is instructive to see the free dispersion relations first. With the free background (i.e. 
Un = 1 everywhere) it is easy to calculate the eigenvalue analytically in momentum space 
to find 



Re(Af ree ) = 1 — 2/cr(cospi + cosp2 + cos»3 + cos p^j , 
Im(Af rcc ) = ±2k^J (sinpi) 2 + (sinp 2 ) 2 + (smp 3 ) 2 + (sin» 4 ) 2 



(4.5) 



Although this expression is valid only for the free background, it turns out to be quite useful 
to understand the eigenvalue distribution in a qualitative level even at strong coupling as 
we will see shortly. 

Usually Re(A) < gives the condition that the Aoki phase appears. In the free case, 
therefore, the Aoki phase has a window |k| > l/(8r), while the Aoki phase condition is 
| | > l/(4r) in the strong coupling limit. Now, as we mentioned in Introduction, it is 
confusing that the diquark condensation has exactly the same criterion for the onset, as 
demonstrated in Figs, [j] and ^. Then, a question arises; which of the diquark superfluid 
phase and the Aoki phase is more favored? The rest of this paper will be devoted to 
answering this question. 




Figure 5: Eigenvalue distribution for a random gauge configuration in the Wilson fermion formal- 
ism at /i = on the 4 lattice for various combinations of k and r. 

Let us see the parameter dependence of the eigenvalue distribution for a random con- 
figuration in the zero density case, which is shown in Fig. |5[ When we increase k with r 
fixed as shown in the lower two figures in Fig. ||, the distribution range is enlarged. We 
can understand this qualitatively from Eq. (|4.5| ) in the free dispersion; both Re(A) — 1 
and Im(A) are proportional to k. The upper two plots in Fig. |B| show the r dependence 
with k fixed. In this case in turn the distribution stretches only along the real axis. This 
feature is also manifest in Eq. ( f4.5| ) since only Re(A) — 1 is multiplied by r. As we can 
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see, the distribution penetrates into the negative real region between (r = 1, k = 0.2) 
and (r = 1, k = 0.3), which is consistent with the known fact that the critical coupling is 
(r = 1, K = 0.25) in the strong coupling limit. In what follows we will employ a value of 
k = 0.23 which is close to the critical point but still outside of the Aoki phase region, if 
any. 

Next, we will investigate how the chemical potential affects the eigenvalue distribution. 
Let us consider the free case first again in which the fourth component is replaced as 
P4 — > p\ — i/i. Then, the real and imaginary parts in the free dispersion are, respectively, 
modified by 



cos(p4 
[sin(p4 — i/i 



i/x) 

2 



cosh fi cos j»4 + i sinh /j sin p^ , 

(sinp4) 2 — (sinh/i) 2 — ^isinh(2/j) sin(2p4) 



(4.6) 



In this simple case it is interesting to see how the free known results are affected by the 
effect of the finite chemical potential which we show in Fig. |(| 
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Figure 6: Eigenvalue distribution for the free Wilson fermion at k = 0.23 on the 4 4 lattice for 
fi = 0, 0.3, 0.6, and 0.9. 

To draw Fig. || we have discretized the momenta pi, p%, P3, and p^ in the range [— 7r, tt] 
into twenty points with equal spacing. We did so in order to make the "density" perceivable 
from Fig. |6| if the momentum is close to continuum with many points, the distribution 
except for the \i = case does not have the empty region strictly. The concentration 
would be hard to see. It is quite interesting to observe a non-trivial structure emerging at 
high n unexpectedly. It is apparent that the density has a similar effect as the hopping 
parameter k; the eigenvalue profile becomes wider in the complex plane as either /i or n 
gets greater. The density modification is not such simple, however, and we presume that 
the rich contents in dense quark matter are attributed in part to this structural difference. 
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At the same time, however, we have to keep in mind that this complicated structure at 
large fi looks like coming from the non-trivial entanglement between different doubler sec- 
tors. In the vicinity of the continuum limit at fi = only the far left edge part corresponds 
to the lightest physical excitation and four other points crossing the real axis are doublers 
going to infinity. This clear separation is missing in view of the eigenvalue distribution of 
Fig. H] at (j, = 0.6 or at \i = 0.9 for instance. This poses a serious question; even though 
we could solve the sign problem somehow, it should be a subtle issue how to separate the 
doublers out at density high enough to allow for excitations of unphysical doublers. 
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Figure 7: Eigenvalue distribution for a random gauge configuration in the Wilson fermion formal- 
ism at k = 0.23 on the 4 4 lattice for = 0, 0.3, 0.6, and 0.9. 



For a randomly generated gauge configuration the fi dependence of the eigenvalue 
distribution reflects the above mentioned structure as displayed in Fig. |7[ Needless to say, 
we can follow the same path to evaluate the chiral condensate but the resulting condensate 
is finite and almost constant independent of the density. This is because naive chiral 
symmetry is explicitly broken in the Wilson fermion formalism, and so we will not present 
the results. Let us now evaluate the diquark condensate in the same way as we did in the 
staggered fermion formalism. It is calculable from the integrated spectral number n£)(£) 
for the operator T>(fi)'D^ (/i). Our results are shown in Fig. ||. 

In the case of the Wilson fermion there are four times degrees of freedom than the 
staggered fermion and so the saturation effect is not yet relevant in the right of Fig. ||. It 
should be mentioned that we measure the slope at the origin in the same way as in the 
staggered fermion; we fit the data up to £ = 0.1 by ud(0 = a £ + b£ 2 . If we change the 
fitting range and the fitting functional form, we would have quantitatively different results. 
The systematic error is not well under control. At least, however, we can state that such 
simple calculations in this work could capture qualitative features of diquark superfluidity. 
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Figure 8: Left) Histogram of no(0- Right) Diquark condensate as a function of p. 



Finally let us discuss the possibility of the Aoki phase from the point of view of the 
spectral density ph(v) corresponding to the parity-flavor breaking condensate. That is 
understood from the 3D plot for the spectral density given in Fig. ||. 

Im(p) ^ 06 Im(p) ^ =09 




Figure 9: Imaginary part of the spectral density, lm(pu (??)) in the complex plane for various values 
of the chemical potential. 

According to the Banks-Casher-type relation we obtained, the parity-flavor breaking 
condensate is to be acquired from the height of lm(pH(il)) at the origin. It should be the 
finite-volume effect that the origin looks smooth and the symmetry breaking looks like 
not occurring even at p = 0. We anticipate that the standing wall would be more sharp 
upright around the origin in the thermodynamic limit. Even in the thermodynamic limit, 
however, the wall has a finite slope at p, ^ which reminds us of the chiral condensate 
discussed in Fig. ^. Thus, the same conclusion can be drawn; the parity-flavor breaking 
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condensate thus takes a non-zero value in the limit of H — * while keeping /i = strictly. 
In the presence of infinitesimal chemical potential, in contrast, the situation changes and 
the condensate is vanishing in the limit of fj, — ► after taking the limit H = 0. Therefore, 
in the exactly same sense as the chiral dynamics we should conclude that there is no 
parity-flavor breaking condensate in this system. As long as the two-color QCD with two 
degenerate flavors is concerned, we do not have to care about the Aoki phase even in the 
Wilson fermion formalism on the lattice. 



5. Remarks 

We saw the eigenvalue distribution of the Dirac operator at finite density T>{fj) and its 
relatives T>(/j,)T)l and V(n)D^ (— //) to discuss the fate of the chiral condensate, the 
diquark condensate, and the parity-flavor breaking condensate. 

We have a conjecture that a similar pattern in the eigenvalue distribution should appear 
also in dense QCD with three colors; the eigenvalue distribution reaches the origin at the 
onset for nuclear matter. At this stage we have no idea what kind of characteristic feature 
is associated with the color superconducting phase. We believe, in principle, that we can 
pursue our strategy in order to access color superconductivity. It has not been successful 
so far to describe the color superconducting phase in the strong coupling limit [51]. Since 



our method does not assume any mean-field nor truncation, our method should be useful 
to clarify what is going on in the diquark channel in strong-coupling QCD. This is what 
we are planning to do as a future extension. 

The present work is focused on the numerical outputs. It is maybe an interesting 
question how the change in the eigenvalue distribution could be interpreted in analogy 
to known phenomena such as the chiral symmetry breaking interpreted as the Anderson 



localization [52, |53|, p% . This research deserves further investigation. Also, it should be 
feasible, in principle, to apply our idea to the weak-coupling regime close to the continuum 
limit using the open gauge configurations if they are available. Since the physical units in 
the color SU(2) world are obscure, unfortunately, the continuum limit is not quite lucid 
then. Nevertheless, in view of the qualitative success of the strong coupling expansion to 



understand hot and dense QCD [48], we may well anticipate that a smooth crossover links 
the strong-coupling regime to the weak-coupling one. This could be checked by inclusion 
of the finite (3 corrections. 

Finally let us mention on the possible extension to the overlap fermion where exact 
chiral symmetry can be defined on the lattice. Then, there is no need to consider the Aoki 
phase from the beginning because the eigenvalue distribution sits on a single circle line 
at [i = 0. This nice feature breaks down, however, at finite density. This is because the 
75-hermiticity is lost at fj, 7^ but it is amazing that chiral symmetry is still realized J38[ j . 
The overlap fermion surpasses the staggered and Wilson fermions; we would be able to 
treat not four but two flavors and look into the behavior of the chiral condensate as well 
as the diquark condensate in the overlap fermion formalism. This extension is also on our 
list for future perspective. 
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